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38 要 : 针对 风云 三 号 C 星 微 波 湿 温 探测 仪 (FY-3CIMMWHTS) 的 陆地 晴空 观测 资料 ， 建 立 了 一 维 变 分 反 演 系统 ， 对 大 
气 的 温 湿 廓 线 进行 反 演 。 为 了 更 好 的 描述 温 湿 廊 线 的 相关 性 , 同时 减 小 温度 和 湿度 在 反 演 过 程 中 相互 之 间 的 误差 传递 ， 
提出 了 使 用 背景 协 方差 矩阵 的 联合 矩阵 和 单独 抢 阵 进行 组 合 反 演 的 方法 。 对 于 MWHTS 模拟 亮 温 和 观测 亮 温 之 间 的 偏 
， 使 用 逐 扫 描 点 的 统计 回归 方法 进行 校正 。 选 择 我 国 部 分 陆地 区 域 的 晴空 观测 亮 温 进行 温 湿 廓 线 的 反 演 ， 并 利用 欧 
H 中 期 天 气 预 报 中 心 (ECMWF) 再 分 析 数 据 、 美 国 国家 环境 预报 中 心 (NCEP) 分 析 数 据 以 及 无 线 电 探 空 观测 (RAOB) 数 据 
对 有 反 演 结果 进行 验证 ， 温 湿 廓 线 的 反 演 结果 与 ECMWEF 再 分 析 数 据 验证 的 最 大 均 方 根 误差 分 别 是 2.59 K 和 11.8796, 

与 NCEP 分 析 数 据 验 证 的 最 大 均 方 根 误差 分 别 是 1.88 K 和 21.50%, 与 RAOB 数据 验证 的 最 大 均 方 根 误差 分 别 是 3.43 K 
和 25.48%， 验 证 结果 表明 了 反 演 结果 的 可 靠 性 。 另 外 与 国外 同类 载荷 AMSU WAE E HAND SAUCE 
度 进行 了 对 比 ， 结 果 表明 : MWHTS RA Be ah i SZ ER 2 UL reg da BE ERZR ERI RE 7J, HEI MWHTS 的 观测 亮 温 
建立 的 一 维 变 分 反 演 系统 具有 较 高 的 反 演 精 度 。 与 NCEP 6 小 时 预报 部 线 的 验证 结果 表明 : 反 演 的 湿度 廓 线 可 以 提高 
预报 廓 线 的 精度 。 
关键 词 : MWHTS， 一 维 变 分 反 演 ， 温 湿 廓 线 ， 逐 扫描 点 校正 ， 预 报 廊 线 
中 图 分 类 号 : TP723,1 文献 标志 码 : A 
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引用 格式 : 贺 秋 瑞 ， 王 振 占 ， 何 杰 颖 . 2016. 用 FY-3C/MWHTS 资料 反 演 陆 地 晴空 大 气温 湿 廓 线 . 


1 引 


的 描述 (Li et al., 2000)。Boukabara 等 人 发 展 的 微波 整 
体 反 演 系 统 (MIRS) 使 用 物理 统计 方法 同时 利用 
AMSU、 微 波 湿度 传感器 (MHS) 以 及 特种 传感器 微波 
成 像 仪 (SSML/S) 的 观测 资料 进行 全 球 大 气 参 数 反 演 ， 


Dll 


温 湿 廓 线 作为 描述 大 气 状 态 的 两 个 重要 参数 ， 在 
E qp 泛 的 应 用 ， % With pco t 
s MS: Put dn peu 该 反 演 系统 具有 全 天 候 的 反 演 能 力 (Boukabara et al., 
(Zavodsky et al., 2012; Araki et al., 2015; Chan et al., UR uu e o o T a AT 
2015). 利用 星 载 迁 感 仪器 探测 大 气温 湿 信息 已 有 超过 。 同 的 先 验 信息 将 会 导致 不 同 的 反 演 精度 。 而 对 于 温度 
50 年 的 历史 ， 传 感 器 的 硬件 设置 取得 了 长 足 的 进步 ， SPEAR UE, Dish EAE Aa RN E TA 
同时 针对 各 种 传感器 的 观测 资料 也 发 展 了 很 多 反 演 算 。 上空 ， 且 在 近 地 层 精度 最 差 ， 其 原因 与 陆地 表面 类 开 
法 。 对 于 微波 波段 的 辐射 计 来 说 ， 王 睫 等 人 和 用 统计 的 复杂 性 以 及 其 上 空 的 温 湿 参 数 变化 特点 有 很 大 关系 
回归 方法 对 先进 微波 探测 装置 (AMSU) 观 测 资料 反 演 EN 
了 海洋 晴空 的 大 气 湿度 廓 线 ， 并 利用 NCEP KIIR OS EM UH 
线 进行 个 例 分 析 ， 一 致 性 较 高 (Wang et aL, 2010). 独立 反 演 温 湿 | 廊 线 的 能 ; 。 目 前 ,关于 FY-3CIMWHTS 
Mathur 等 人 同样 利用 统计 回归 方法 对 微波 湿度 探测 。 ”的 硬件 设计 及 数据 定 标 验证 已 有 相关 文献 进行 介绍 
SAPHIRA ES Ls A o A db eh i 
d Ru uM et al., pod Karbou 反 演 是 通 ; — AR HERE BOR 射 传输 方 
等 人 使 用 神经 网 络 的 方法 对 AMSU 观测 资料 进行 陆 。“” 程 求 逆 米 实现 的 ， 是 提高 反 演 精度 的 根本 途径 。 因 此 
地 晴空 大 气温 湿 廓 线 反 演 ， 针 对 陆地 表面 类 型 和 大 气 ”本文 使 用 物理 有 反 演 方法 对 MWHTS 陆地 晴空 观测 完 温 
参数 的 季节 变 化 进 行 了 分 类 处 理 ， 并 得 到 较 高 的 精 度 进行 反 演 看 JL o 然而 ， 由 于 在 大 气 辐射 传输 建 模 的 过 
(Karbou et al., 2005). Gangwar 等 人 使 用 同样 的 方法 对 程 中 ， 陆 地 表面 发 射 率 的 计算 难度 大 且 计 算 精 度 差 


^F 


AMSU 观测 资料 进行 海洋 和 陆地 的 晴空 EA HR EEA Fe al., VENT i et al, A 所 以 本 文 为 


cA pe ucc mE 详细 分 类 外 BE, 而 是 从 先 验 信息 的 有 度 出 发 来 提高 反 
廓 线 分 类 反 演 ， 并 对 反 演 误差 进行 了 详细 的 分 析 PORE d 
(Gangwar et al, 2014). Li 等 人 发 展 的 先进 垂直 探测 器 i LM ME cleri rn ann sie 


数据 处 理 软件 包 (IAPP) 对 AMSU 和 高 分 辩 率 红外 辐射 再 先 代 的 一 维 变 分 肥 演 系统 ， 对 中 国 部 分 陆地 区 域 旱 
探测 仪 (HIRS) 观 测 资 料 进行 了 有 云 和 无 云 条 件 下 大 气 空 条 件 下 的 大 气温 湿 廓 线 进行 了 反 演 研究 。 利用 2015 


; —— 年 4 月 和 5 月 两 个 月 的 亮 温 数据 进行 反 演 ， 并 分 别 与 
参 ER, Ex [y Jl LEE E zx EIE J 2 
RR A E E E RUM ECMWF 再 分 析 数 据 、NCEP 分 析 数 据 以 及 RAOB 数 
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据 进行 了 对 比 验 证 分 析 , 同时 也 与 同类 载荷 AMSU 观 
测 资料 的 反 演 精度 进行 了 对 比 。 


2 资料 与 模式 


2.1 MWHTS 资料 


MWHTS 是 超 外 差 式 接收 机 的 全 功率 型 微波 辐射 
计 ，8 个 大 气温 度 探测 通道 设置 在 118.75 GHz 氧气 吸 
收 线 附 近 ，5 个 大 气 湿度 探测 通道 设置 在 183.31 GHz 
水 汽 吸收 线 附近 ， 两 个 窗 区 通道 分 别 设置 在 89 GHz 
和 150 GHz。 设置 在 118.75 GHz 的 这 组 探测 通道 是 国 
际 上 业务 卫星 首次 使 用 的 大 气 探测 通道 ， 与 183.31 
GHz 通道 联合 可 以 实现 大 气温 湿 参 数 的 同时 探测 。 
MWHTS 扫描 宽度 约 2700 Km， 星 下 点 空间 分 辨 率 约 
15 Km， 每 条 扫描 线 有 98 个 扫描 点 ， 每 个 扫描 点 均 对 
光一 个 卫星 天 顶 角 ， 扫 描 周 期 2.667 秒 。 另 外 ， 关 于 
MWHTS 的 通道 权重 函数 随 高 度 的 分 布 特性 ，Bao Xt 
行 了 详细 的 分 析 (Bao, 2014). 


2.2 研究 资料 与 模式 


本 研究 中 所 使 用 的 数据 有 : (1) FY-3C/MWHTS 
亮 温 数据 ， 由 中 国 气 象 局 网 站 L1 数据 产品 提供 ， 时 
间 为 2014 年 2 月 到 2015 年 5 A; (2) NOAA-18 星 
AMSU-A 和 AMSU-B 亮 温 数据 ， 由 中 国 气象 局 网 站 
Llc 数据 产品 提供 ， 时 间 为 2014 年 2 月 到 2015 年 $ 
H; (3) ECMWF 再 分 析 数 据 (ERA Interim)， 包 括 温 湿 
亡 线 以 及 表面 参数 。 其 中 温 湿 廊 线 从 地 面 到 高 空 共 分 
为 37 JE, 包括 1000, 975, 950, 925, 900, 875, 850, 
825, 800, 775, 750, 700, 650, 600, 550, 500, 450, 
400, 350, 300, 250, 225, 200, 175, 150, 125, 100, 


70, 50, 30, 20, 10, 7, 5, 3, 2 和 1 hPa， 分 辩 率 
是 | 0257x025 , 由 ECMWF W 站 
(http://apps.ecmwf.int/datasets) 提 供 ， 时 间 为 2014 年 2 
月 到 2015 年 5 月 ; (4) NCEP 全 球 预报 系统 (GFS) 6 
小 时 预报 数据 和 分 析 数 据 , 包括 温 湿 廓 线 和 表面 参数 ， 
其 中 温 湿 廓 线 从 地 面 到 高 空 共 分 为 26 层 , 包括 1000, 


975, 950, 925, 900, 850, 800, 750, 700, 650, 600, 
550, 500, 450, 400, 350, 300, 250, 200, 150, 100, 
70, 50, 30, 20 和 10 hPa, 分辨 率 是 0.25°x0.25°, H 
NCEP 网 站 (http://rda.ucar.edu/datasets/ds084.1) 提 供 ， 
时 间 为 2015 年 1 月 到 2015 ££ 5 A, XT- ECMWF 
再 分 析 数 据 的 廓 线 分 层 相 一 致 ， 对 NCEP 廓 线 进行 三 
次 样 条 内 插 处 理 ; (5) RAOB 数据 ， RAOB 网 站 
(http://www.esrl.noaa.gov/raobs) Heft, 时 间 为 2015 年 4 
月 到 2015 年 5 A, 同时 对 RAOB 数据 进行 质量 控制 ， 
温 湿 探 测 总 的 压强 层 数 大 于 20 且 最 小 压强 值 大 于 200 
hPa 为 有 效 值 ， 同 样 使 用 插值 的 方法 处 理 成 37 层 。 因 
为 本 文 是 针对 中 国 部 分 陆地 区 域 上 空 大 气 参数 的 反 演 
研究 ， 所 以 以 上 数据 的 地 理 范 围 均 选择 为 
(25°N~45°N,90°E~120°E)。 对 于 模拟 亮 温 的 计算 ， 本 
文 使 用 的 是 ECMWF 发 展 的 快速 辐射 传输 模型 
(RTTOV v11.2) (Hocking et al, 2014). 


2.3 晴空 数据 的 选择 


关于 晴空 数据 的 选择 ， 己 有 很 多 文献 进行 了 研 
究 。Karstens 等 人 通过 设置 相对 湿度 的 门限 值 来 进行 
晴空 的 判断 ,但 是 水 汽 参数 随时 间 和 空间 的 变化 较 大 ， 
因此 该 门限 值 有 一 定 的 主观 性 ， 在 具体 应 用 时 需要 根 
据 当 地 气象 条 件 进 行 统计 分 析 得 到 (Karstens et al., 
1994)。Ishimoto 等 人 利用 红外 云 产 品 与 微波 亮 温 数据 
进行 匹配 来 进行 晴空 数据 的 选择 ， 但 是 该 方法 对 红外 
云 产 品 的 依赖 性 很 大 (Ishimoto et al., 2014). Buehler 
等 人 利用 微波 传感器 各 个 通道 之 间 的 亮 温 数据 关系 进 
行 云 雨 数据 的 选择 ， 但 是 对 薄 云 区 域 的 亮 温 数据 不 能 
iQ oU QU t Ret ee. 
确 的 描述 晴空 大 气 ， 避 免 云 带 来 发 射 、 吸 收 以 及 散射 
的 影响 ， 根 据 ECMWF 再 分 析 数 据 中 的 云 水 路 径 为 0 
来 选择 晴空 大 气 参数 ， 然 后 与 MWHTS 和 AMSU 亮 
温 数 据 进行 时 间 和 空间 上 的 匹配 , 选择 晴空 观测 数据 。 


3 一 维 变 分 反 演 系统 


3.1 一 维 变 分 反 演 算法 原理 


唱 射 传输 方程 是 一 个 非 线 性 方程 , 大 气 参数 通过 
对 其 求 逆 得 到 ， 因 此 这 又 是 一 个 欠 定 的 求解 过 程 。 一 
维 变 分 反 演 算法 是 物理 反 演 方法 中 的 典型 代表 ， 该 算 
法 主要 包括 两 部 分 : 一 个 是 计算 正 演 亮 温和 亮 温 梯度 
的 辐射 传输 模型 ， 另 一 个 是 用 来 权衡 先 验 信息 和 卫星 
观测 亮 温 对 最 终 解 页 献 的 最 小 化 代价 函数 。 假 设 观 测 
误差 和 先 验 信息 误差 是 无 偏 和 不 相关 的 ， 并 且 服 从 高 


斯 分 布 ， 那 么 大 气 参 数 x 的 最 优 估计 可 以 通过 最 小 化 
代价 函数 (Boukabara et al., 2011): 
J = (2-8 B œ- +H) - (1) 


I| R"[H(9) - I] 
来 求解 。 其 中 是 背景 廊 线 ，B 是 背景 协 方差 矩阵 ， 
R 是 测量 误差 协 方差 矩阵 ， 包 括 观 测 误差 协 方差 矩阵 
E 和 正 演 误差 协 方差 矩阵 F, I 表示 观测 亮 温 ，H(x) 
表示 对 大 气 参 数 x 计算 模拟 亮 温 , T 表示 和 矩阵 的 转 置 。 
对 代价 函数 求 导 可 得 ; 
VJ -B'(x-x)-H'QO)R"[H(x)- I] (2) 
其 中 ,五 为 正 演 亮 温 的 梯度 ， 它 表示 正 演 亮 温 值 对 大 
SB x 值 的 变化 的 灵敏 度 。 通 过 让 式 子 (2) 为 零 来 求 
得 最 优 解 ， 使 用 牛顿 迭代 算法 可 得 : 
x, =x’ + BH" (x JH (x )BH (x,)+ 
R|'U - H(x)- H (x, Xx" —x,)] 
其 中 ,nn RIERO xi 表示 式 子 (2) 的 最 优 解 ， 
即 反 演 廓 线 。 从 式 子 (3) 中 可 以 看 出 先 验 信 息 B x? 
以 及 观测 亮 温 和 模拟 亮 温 之 间 的 偏差 TH(x,) 对 反 演 精 
度 有 直接 影响 。 


3.2 先 验 信息 


求解 欠 定 的 非 线 性 方程 , 需要 在 无 数 个 解 中 寻找 
最 优 解 ， 且 该 解 符合 实际 的 物理 意义 ， 必 须 使 用 先 验 
信息 进行 限制 ， 它 直接 决定 了 解 的 存在 性 和 精度 。 背 
景 协 方差 矩阵 如 使 用 气候 学 数据 集 产生 , 包括 RAOB 
数据 , 同化 系统 的 再 分 析 数 据 等 (Rodgers, 2000)。 本 文 
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使 用 2014 年 一 年 内 ECMWF 再 分 析 数 据 : 


下 的 111069 条 温 湿 亡 线 生 成 B， 具 体 公式 为 
1 


oj = 33e xo, 


1 j=l 


25 x Ay FN AS AW hE EB Be AL HA PE JR AG IY 
Ti d BBR HZ 
单独 矩阵 ， 当 x* 同时 
成 的 是 可 以 描述 温 湿 廓 线 相互 间 的 相关 性 
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(a) Individual matrix 


温 湿 廓 线 分 层 数 
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温 湿 廓 线 分 层 数 

(b) 联合 矩阵 
(b) United matrix 
1 相关 系数 矩阵 


Fig. 1 Correlation matrix 
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晴空 条 件 


(4) 


其中 ,oy 表示 背景 协 方差 矩阵 中 第 i 行 第 j 列 的 元 素 ， 
x 表示 廊 线 x 的 均值 ，N 是 所 使 用 的 廓 线 样本 数目 。 
生成 的 是 
间 的 相关 性 的 协 方差 矩阵 ， 称 为 
示 温 度 廓 线 和 湿度 廓 线 时 ， 生 
的 协 方差 矩 
E 阵 更 
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其 中 ， 温 湿 廓 线 分 层 数目 1-37 表示 从 高 空 到 地 


外 的 温度 廊 线 的 相关 系数 ，38~74 表示 从 高 空 到 地 面 
的 湿度 廊 线 的 相关 系数 。 从 图 1(a) 可 以 看 出 ， 温 度 和 


湿度 在 各 自 层 之 间 具 有 很 强 的 相关 性 ， 但 互相 之 间 不 
存在 相关 性 。 而 从 图 1(b) 可 以 发 现 ， 温 度 和 湿度 除了 


各 自 层 之 间 的 强 相关 性 外 ， 它 们 之 间 互 相 性 很 强 ， 特 


别 是 在 高 空 时 ， 互 相关 系数 达到 了 0.8 以 上 。 在 自然 


条 件 下 ， 
时 它们 对 相互 之 间 的 探测 通道 也 有 辐射 贡献 。 


大 气 的 自然 变化 状态 ， 有 助 于 提高 反 演 精度 。 


大 


用 联合 矩阵 对 温 湿 廓 线 进 行 反 演 才能 更 加 准确 


日 
在 反 演 过 程 中 这 两 个 参数 的 不 确定 性 将 会 带 来 相 
响 ， 特 别 是 对 湿度 反 演 来 说 ， 如 果 不 能 够 提供 
确 的 先 验 温 度 廓 线 ， 那 么 这 将 使 湿度 反 演 的 精度 大 大 


温度 廓 线 和 湿度 廓 线 本 身 就 存在 相关 性 ， 同 
比 使 


述 


互 影 
足够 精 


减 小 (Aires et al., 2015)。 在 前 人 建立 的 一 维 变 分 反 演 系 


都 是 使 用 背景 协 方差 矩阵 的 单独 矩阵 或 联合 矩 
并 未 有 人 同时 使 
., 2011; 
陆地 上 


二 者 进行 反 演 研究 (Boukabara 
本 文 针 对 MWHTS 的 


Li Jun et al., 2000). 


a Sci, Diui SEAR AE EA 
进行 组 合 反 演 的 方法 ， 即 反 演 温度 廓 线 时 ， 使 
矩阵， 舍弃 此 次 反 演 得 到 的 湿度 亡 线 ， 而 反 演 湿度 廓 
线 时 ， 使 ) 
次 反 演 分 别 得 到 了 大 气 的 温 湿 
充分 不 


单独 矩阵 
联合 


单独 矩阵， 舍弃 此 次 反 演 的 温度 廊 线 ， 两 
RAB, PEA TEBE 
TERRAE, E X REN SCIBUSE 


REP eee fe ie RA » 


背景 
HIR 


JT ERA x? 来 说 ， 
亡 线 对 反 演 结果 有 直接 影 
越 有 利于 提 
气候 学 数据 集 的 平均 廓 线 ; 
(3) 时 间 和 空 
平均 廊 线 只 是 对 该 区 域 历史 探测 数据 求 均值 ， 


| 式 子 (3) 中 可 以 看 出 ， 
向 , 其 越 接近 真实 廓 线 ， 
高 反 演 精度 。 双 可 使 用 的 数据 源 包括 : (1) 
(2) 统 计 反 演 的 输出 廓 线 ; 
3 间 上 匹配 的 预报 廊 线 。 因 为 气候 数据 集 的 
与 真实 


大 气 状态 偏差 较 大 ， 而 当 使 用 统计 反 演 方法 的 输出 廓 


线 时， 需要 
计算 ， 


明显 增加 计算 量 。 


首先 进行 统计 样本 的 选择 以 及 回归 系数 的 
因此 本 文选 取 NCEP 6 小 时 


预报 廓 线 作 为 背景 廓 线 。 同 时 把 反 演 结果 与 预报 廓 线 


的 精度 进行 对 比 , 来 说 明 MWHTS 观测 资 
气 预 报 同 化 系统 的 意义 。 而 对 于 初始 猿 想 廊 线 ， 
子 (3) 中 xo 的 值 来 说 , 7 


资料 对 数值 天 
即 式 
需要 从 概念 上 与 背景 廓 线 进行 区 


分 ,一 个 好 的 初始 猜想 廊 线 可 以 加 速 迭代 过 程 的 收敛 ， 


但 对 反 演 精度 3 
统 中 以 


无 影响 ， 本 文 构建 的 一 维 变 分 反 演 系 


背景 廓 线 代 蔡 初始 猜想 廓 线 。 


3.3 FY-3C/MWHTS 通道 偏差 校正 


与 海洋 表 


看 发 射 率 的 计算 相 比 , 陆地 表面 发 射 率 


的 计算 误差 较 大 ， 从 而 导致 了 不 同 于 海洋 情况 下 模拟 


亮 温 的 系统 偏差 。 该 系统 偏差 直接 影响 反 演 精度 。 本 
文选 择 2014 年 2 月 到 2015 *E 3 H 14 个 月 的 晴空 条 件 
下 ECMWF 再 分 析 数 据 进行 模拟 亮 温 的 计算 ， 其 地 


表 发 射 率 使 用 快速 发 射 率 模 型 版 本 5(FASTEM-5) 进 行 
计算 ,为 了 对 MWHTS 所 有 扫描 点 的 统计 特性 进行 分 


pr, B 


分析 数据 与 MWHTS 观测 亮 温 的 匹配 规则 设置 
EE OU 
共 得 到 了 超过 十 万 个 匹配 对 。 


经 纬度 误差 分 别 小 于 0.05*, 
通过 对 模拟 亮 温 和 观测 


^c lbi 


进行 统计 分 析 发 现 大 部 分 通道 具有 线性 统计 关 


系 , 计算 每 个 通道 中 模拟 亮 温和 观测 亮 温 的 相关 系数 ， 
其 随 扫描 点 的 变化 如 图 2 所 示 。 


表 发 射 率 计 算 误差 的 影响 与 观测 亮 温 2 
大 , 其 大 部 分 扫描 点 的 相关 系数 较 低 ， 
通道 2 的 相关 系数 在 所 有 扫描 点 处 均 较 低 ， 保 持 在 
0.75 以 下 ， 这 可 
度 廓 线 精度 不 够 高 有 关 。 


图 2 可 知 ， 


窗 区 通道 1 和 10 的 模拟 亮 温 受 地 


间 的 误差 较 
都 在 0.85 以 下 。 


能 与 在 高 空 ECMWF 再 分 析 数 据 的 温 
对 于 其 余 温 度 通道 来 说 ， 通 


道 4、5、6、7 具有 较 高 的 相关 系数 ， 均 保持 在 0.95 


附近 ， 


而 通道 3、8、9 相关 系数 稍 低 ， 但 大 部 分 扫描 


点 处 均 保持 在 0.90 以 上 。 湿度 通道 的 相关 系数 均 高 于 


0.80， 


if RIA 


了 较 大 的 变化 ， 大 的 卫星 天 项 角 


时 均 保 持 在 0.90 附近 ,而 小 的 卫星 天 顶 角 时 保持 在 0.8 


附近 。 


总 体 来 说 ， 每 个 通道 的 相关 系数 都 随 扫描 点 有 


一 定 的 起 伏 ， 除 了 数据 时 空 匹 配 带 来 的 误差 外 ， 这 是 


€ 


才 


j 交 轨 扫 描 方式 的 传感器 所 获得 观测 数据 的 特点 ， 


该 扫描 方式 获得 的 观测 数据 往往 对 卫星 天 顶 角 有 具有 依 和 矩阵。 利用 各 个 通道 校正 后 的 亮 温 与 模拟 亮 温 之 间 的 
赖 性 。MWHTS 十 五 个 通道 的 模拟 亮 温 和 观测 亮 温 的 偏差 求 平方 得 到 正 演 误差 协 方差 矩阵 E 的 对 角 元 素 ， 
平均 偏差 随 卫 星 天 顶 角 的 分 布 如 图 3 所 示 。 j 对 于 观测 误差 协 方差 矩阵 下 的 对 角 元 素来 说 ,使 用 
E UR _ A ut 通道 的 灵敏 度 的 平方 值 来 产生 dshimoto et al., 2014), 
Ss 本 文 使 用 的 是 Bao X} FY-3C/MWHTS 在 轨 运 行 数据 所 
计算 灵敏 度 的 平均 值 (Bao，2014)。 其 中 正 演 误差 和 通 
道 灵 敏 度 的 取 值 如 表 1 所 示 。 
#1 正 演 误差 与 通道 灵敏 度 


br 
0.9 ] S7 
5 


0.8 ] 
i: Table 1 The forward model error and channel sensitivity 
& 通道 IENOROUK ARREK 通道 正 演 误差/K 通道 灵敏 度 /K 
0.7 1 8.30 0.23 9 2.36 0.27 
2 2.39 1.62 10 4.02 0.34 
3 0.93 0.74 11 1.50 0.47 
0.64 4 0.82 0.59 12 1.06 0.34 
5 0.85 0.65 13 2.89 0.30 
! T T T | 6 1.06 0.52 14 1.66 0.22 
" a ur id 100 7 1.22 0.49 15 1.54 0.27 
DERI 首 2 通道 3; MÉ 4 MĚ s 8 2.12 0.27 
Ll Hio; -- -通道 7; -通道 8 ----- Eo; -通道 10 m 
-- 通 道 11; 38b 12 - -通道 13; ---- 通道 14， 通道 15 3.4 质量 控制 
图 2 观测 亮 温 和 模拟 亮 温 的 相关 系数 x 
Ml H 
Fig. 2 The correlation between observed and simulated AS SOIT i IE (OS Stbs iE AML MUS 亮 温 质量 控 
rightness temperature 制 标准 来 对 反 演 结果 进行 质量 控制 。 选 择 的 迭代 收敛 
i ] i 标准 是 当代 价 函 数值 的 相对 变化 范围 在 0.01 PN TEA 
a oh ET irr "Pa pes 停止 ， 且 最 大 迭代 次 数 为 10， 由 于 超过 10 Ses 
z8 zB re B XA 
W $05 50 100 28 8 50 100 #5 ^ 50 100 AE, 背景 廊 线 BUI Ay Bz yes BR ZR o 利用 背景 廓 线 的 模 
扫描 点 扫描 点 Hd AT 
€ oa ££ poo 拟 亮 温 与 观测 亮 温 的 差 让 对 输入 的 背 景 廓 线 和 观测 亮 
moy UU a? 通道 6 JH 3 
& 4 m wk CERE Ge 50 100 温 进 行 质量 控制 ， 如 果 任 一 通道 的 差 值 大 于 20K, 8 
H1 iili 点 m aA m EE 弃 该 组 亮 温 数据 。 
= 5[F wer] x $e 关 10 [一 通通 可 
$ o = NM 8p * 
2 56 50 100 2 50 50 100 A % 50 100 
m RUM ^ SEI VAS N A. 
ER n pe 4 i iN BBG TAS E 
= oh SEO] & TTE E 
ak : an i^ 19 ， 100 * ^o : sA 100 41 反 演 步骤 
i er aS a E xc ae 通通 后 
Hol cu "PEE amua | PES 本 文 建立 的 一 维 变 分 反 演 系统 反 演 的 基本 步骤 
w “o 50 100:8 ^0 50 — 109 0 50 100 包括 
扫描 点 扫描 点 扫描 点 : 
图 3 MWHTS 通道 亮 温 偏差 对 扫描 角 的 依赖 性 1) 对 FY-3C/MWHTS 进行 晴空 条 件 下 亮 温 数据 
Fig. 3 Scan-angle dependence of MWHTS brightness temperatures 的 选择 ， Ail} j 逐 扫描 点 统计 回 归 的 校正 方 法 并 进行 偏 


bias 
从 图 3 中 可 以 看 出 , 每 个 通道 平均 亮 温 偏差 对 卫 ZE. ee |" 
ERARE- EMRE, HRB KBE, PA 2) 用 校正 过 的 亮 温 数据 与 ECMWF 再 分 析 数 
了 较 大 的 偏差 且 随 天 顶 角 变化 很 大 。 对 于 温度 通道 据 、NCEP 预报 数据 和 分 析 数 据 、RAOB 数据 进行 空 
说 ， 通 道 2、3、4、5、6 的 平均 亮 温 偏 差 随 天 项 角 起 。 ”同和 时 间 上 的 匹配 。 其 中 ，MWHTS 亮 温 数据 与 
伏 不 大 ， 但 是 通道 7、 dd 通道 起 伏 趋势 明显 。 ECMWF 再 分 析 数 据 、NCEP 预报 数据 以 及 分 析 数 据 
对 于 湿度 通道 来 说 , 通道 11、12、13 平均 亮 温 偏差 随 的 匹配 规则 均 采 用 时 间 误 差 小 于 0.5 小 时 ， 经 纬度 误 
sl E Du di 差分 别 小 于 0.05"。 而 与 RAOB 数据 的 匹配 规则 是 时 
明显 。 在 物理 反 演 系统 中 ， 对 于 系统 偏差 进行 校正 万 间 误 差 小 于 3 小 时 ， 经 纬度 误差 分 别 小 于 1°. 
为 重要 ， 否 则 可 能 会 导致 迭代 的 不 收敛 或 者 是 得 到 错 3) 把 匹配 到 的 NCEP 6 小 时 预报 廓 线 作 为 背景 
误 的 反 演 结果 。 本 文 从 统计 学 的 角度 提出 了 逐 扫 描 点 “” 廓 线 、 第 3 节 建 立 的 测量 误差 协 方差 抢 阵 及 和 背景 协 
统计 回归 校正 方法 对 MWHTS 各 个 通道 的 观测 亮 温 进 方 郑 惩 阵 妃 的 单独 惩 阵 和 联合 矩阵 输入 到 一 维 变 分 反 


行 偏差 校正 。 为 了 充分 利用 MWHTS 观测 数据 ， 对 模 PRR. 


POP 


拟 亮 温和 观测 亮 温 的 相关 系数 大 于 0.80 的 扫描 点 进行 D AABE Er A ATE RUIT TR 
偏差 校正 ， 表达 式 为 (Li et al., 2000): MARG CARS AR 合 和 矩阵 进行 反 演 计算 ， 取 使 联合 
I aT, +h, (5) EVP RE ERAS FE ee scii DIRE HG 

M "i a 线 ， 得 到 反 演 结果 ， 并 利用 匹配 到 的 ECMWF 再 分 析 

bb, THER, T ARBETE, a ERE 数据 、NCEP 分 析 数 据 和 RAOB 数据 对 反 演 结果 进行 


b RUE. JUP i RR MWHTS 的 通道 ， 取 值 范围 为 精度 疮 证 分 析 4 

"uo he R 本 文选 择 2015 年 4 月 到 2015 年 5 月 两 个 月 的 
HI EDERE MWHTS 各 通道 之 间 相 关 性 的 情 。 FY-3C/MWHTS 陆地 晴空 亮 温 数据 进行 反 演 计算 。 

况 下 可 以 使 用 和 矩 阵 R 的 对 角 元 素 作为 测量 误差 协 方差 


2 逐 扫描 点 统计 回归 偏差 校正 效果 et o s 


yE pF Haya 
FY-3C/MWHTS 晴空 亮 温 与 ECMWF BAIN "s. IU i Ris. uta DLE 
时 间 误 差 小 于 2 小 时 , 经 纬度 误差 分 别 小 于 0.05? BUG iic) PM 
配 规则 ， 一 共 得 到 10078 个 匹配 对 。 对 匹配 对 中 的 亮 ‘el ug | oda S 
EL MH 一 fea + £ 
温 进 行 偏差 校正 ， 如 图 4 所 示 。 NS Ki 
2005 + 4 200 E. 
A 0 d x A V rir uw x OEA Ae T 300 - s 二、 J 300 b x 
we P05 50 10028 -60 0 1008 70 50 100 a T 
ES 扫描 点 E 44907 i 下 * 
E : 通道 4 £i VM Es 通道 6 EM 5 Eu i 
= op M 通道 5 "M S + x * 
# So 50 100 2 to 50 10049 3g 50 100 = eoo. x = 600 上 i 
" Hi js u 扫描 点 M 习 描 点 Y xk 
= ; ci < ; r 关 700 | + 700} * 
E So E i 100 z 5o 3 To% ES i 100 nii E | ul b 
TES 扫描 点 描 点 F A 
& 20 a Z2 —r g5 à R^ 900 + E 900 + EN E 
E 通道 10 = 通关 11 = 通道 12 t EN 
i: 0 i 0 í Wo $ Ng BW 
w 50 100: 20 50 10058 50 50 100 PM 95 "05 1 12000 1 2 3 ^4 
Hii ah 扫描 点 Hi "T ŽNK) 温度 均 方 根 误差 /(K) 
= ] aur F Es ELE + AAI pen 联合 矩阵 
E Pa E 通道 14 | EE d T (a) 温度 
we 20 50 10038 9g 50 100 8 °o 50 100 (a) Temperature 
Hd ux 扫描 点 SEE — T c - 
一 一 : 校正 前 的 平均 亮 温差 ; ——: 校正 后 的 平均 亮 温差 300- $ J 300 L A 
图 4 通道 的 亮 温 偏差 校正 结果 EY € 
Fig.4 The correction of brightness temperatures bias 400 Ne 4 400 + 
从 图 4 可 以 看 出 窗 区 通道 1 和 10 只 有 部 分 扫 SS 
点 进行 了 偏差 校正 ， 通 道 2 由 于 相关 系数 低 而 未 进行 T JR N 
偏差 校正 。 对 于 其 余 温度 通道 来 说 ， 校 正 效果 比较 明 F om) L | fw " 
显 ， 但 是 通道 7、8、9 在 个 别 扫描 点 处 仍然 存在 较 大 z + 5 E 
的 偏差 。 对 于 湿度 通道 来 说 ， 通 道 12、14、15 校正 效 dini F 7: il 
果 较 好 ， 通 道 11 在 天 底 点 处 校正 效果 不 好 ， 通 道 13 M Hl & 
在 大 的 天 顶 角 处 的 校正 效果 较 好 ， 而 在 小 的 天 顶 角 处 X 
ERN i se aei n 3. * * 
校正 效果 不 理想 。 总 的 来 说 ， 逐 扫描 点 统计 回归 校正 900 上 "p 900 | + 
方法 的 校正 效果 受 通道 模拟 亮 温和 观测 亮 温 之 间 的 相 | * 
关系 数 的 影响 , 因此 在 实际 应 用 中 也 有 一 定 的 局 限 性 。 15 40 5 0 5 5 w s 20 
湿度 偏差 /(%) 湿度 均 方 根 误差 /(%) 
4.3 反 演 结果 的 验证 分 析 (AH. i BAER 
(b) 湿度 
根据 FY-3C/MWHTS 通道 权重 函数 特性 分 析 (b) pr 
Me s 背景 al 
(Bao, 2014), 选取 温度 反 演 廓 线 的 压强 范围 是 10~1000 图 5 背景 协 方差 矩阵 对 反 演 结果 的 影响 
Fig.5 The effect of background covariance matrix on inversion 


hPa。 温 度 反 演 廊 线 的 压强 范围 是 250~1000 hPa. AS 


results 


文 用 平均 偏差 和 均 方 根 误差 作为 对 反 演 结果 统计 验证 432 反 演 结果 的 验证 分 析 
IRE rss (Mathur etat. 2015). 2 利用 地 理 坐 标点 (90.54*E,43.03*N) 的 2015 年 4 
4.3.1 背景 协 方差 秆 阵 对 反 演 结果 的 影响 H 26 号 05:45 UTC 亮 温 单 点 反 演 结果 和 地 理 坐 标点 


分 别 使 用 Hed 方差 矩阵 的 联合 矩阵 和 单独 阜 阵 — (90 sQ»E, 43.00°N) IY 2015 年 4 月 26 号 06:00 UTC 的 


组 亮 温 与 背景 — 亮 温 的 偏差 超过 20 K 而 舍弃 ， 从 图 6 可 以 看 出 ， 反 演 的 温 湿 廊 线 与 ECMWE 

得 到 9483 组 温 湿 廊 线 。 并 利用 ECMWF 再 分 析 数 据 。” 温 温 廓 线 在 结构 上 一 致 性 较 高 ，ECMWTF 温 湿 廓 线 结 

对 反 演 结果 进行 统计 特性 验证 分 析 ， 如 图 5 所 示 。 构 上 小 的 变化 ， 反 演 廊 线 也 可 以 进行 很 好 的 描述 。 同 
从 图 5) 中 可 以 看 出 , 无 论 使 用 单独 矩阵 还 是 联合 。 时 可 以 发 现 温 湿 廓 线 的 反 演 偏差 和 背景 偏差 在 结构 上 


息 阵 ,温度 反 演 廓 线 的 二 沟 俩 差 均 保持 在 09 玉 之 和 内， 随 压 强 有 相同 的 变化 ， 因 此 背景 偏差 可 能 是 反 演 误差 
均 方 根 误差 在 接近 地 面 处 最 大 , 达到 3.9 K。 使 用 联合 的 一 个 来 源 。 

和 矩阵 的 反 演 结果 在 300~1000 hPa 范围 内 , 精度 高 于 使 另外 ， 时 间 匹 配 误差 直接 影响 反 演 精度 的 验证 ， 
用 单独 怎 阵 的 反 演 结果 ， 但 是 在 200-300 hPa 范围 ， 因此 把 时 间 匹 配 误差 缩小 为 0.5 小 时 ， 得 到 105 AA 
使 用 单独 矩阵 高 于 联合 矩阵 。 从 图 5(b) 中 可 以 看 出 使 。 温 数 据 ， 使 用 ECMWF 再 分 析 数 据 对 反 演 结果 进行 验 
j 不 同 抢 阵 的 湿度 反 演 廓 线 的 平均 偏差 基本 上 保持 在 证 分 析 。 利 用 0.5 小 时 匹配 误差 的 未 校正 亮 温和 校正 
OWA, SITAR 7AA, ERRAR 亮 温 进 行 反 演 ， 其 中 分 别 有 9 组 和 1 组 的 观测 亮 温和 
E 阵 进行 反 演 的 精度 在 各 个 压强 处 均 高 于 联合 第 阵 ， 背景 廓 线 的 模拟 亮 温 的 偏差 超过 了 20 K， 因 此 分 别 得 
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Fig. 7 The mean dif 


从 图 7(a)! 


差 相 比 校正 前 有 了 很 大 的 减 小 ， 
E nh 


有 0.1 K 的 增 大 ， 


图 7 反 演 廊 线 的 平均 偏差 和 均 方 根 误差 


ference and root mean square difference of 


inversion profiles 
可 以 看 出 , 偏差 校正 后 的 温度 平均 偏 
最 大 处 可 减 小 2 Ke 
相 比 校正 前 约 
与 温度 探测 通道 7、8、9 的 校正 


效果 有 直接 的 关系 。 在 1000 hPa 附近 校正 前 后 的 反 注 
精度 相当 ， 这 很 大 程度 上 是 因为 窗 区 通道 的 校正 效果 


不 理想 


比 校正 前 均 有 提高 ,其 中 
从 图 7(b) 中 可 以 看 出 ， 
正 前 均 有 减 小 ， 虽 然 通 道 11 IS 
理想 ， 但 是 均 方 根 误 差 在 湿度 探测 通道 起 主要 贡献 作 
用 的 400~750 hPa 范围 


导致 的 ， 而 在 其 他 压强 处 校正 后 的 反 演 精度 相 


， 最 大 均 方 根 误差 为 2.59 Ko 
偏差 校正 后 的 平均 偏差 相 比 校 
偏差 校正 效果 不 


cr 


明显 的 减 小 。 而 在 250-400 


hPa 范围 , 偏差 校正 对 湿度 反 演 精度 基本 上 没有 影响 ， 
其 中 ， 最 大 均 方 根 误差 为 11.87%。 


通过 单 点 反 演 结果 的 对 比分 析 以 及 104 组 反 演 
结果 的 统计 验证 分 析 可 以 发 现 ， 背 景 廓 线 的 误差 与 反 
演 廓 线 的 误差 总 体 趋 势 一 致 性 较 高 ， 因 此 可 以 判断 背 


景 廓 线 的 精度 对 反 演 精度 有 很 大 的 影响 。 
4.3.3 Ej NCEP 分 析 数 据 和 RAOB 数据 的 验证 分 析 


为 了 保证 反 演 


结果 的 可 靠 性 ， 分 别 与 NCEP 分 析 


数据 和 RAOB 数据 进行 精度 验证 分 析 ， 根 据 4.1 节 中 


反 演 步骤 的 匹配 规则 , 分 别 得 到 了 104 和 1068 组 反 演 
结果 ， 验证 结果 如 图 8 所 示 。 
从 图 8 可 以 看 出 ， 对 于 NCEP 分 析 数 据 的 验证 结 


果 来 说 ， 温度 2 


Hn ECMWF 再 分 析 数 据 的 验证 


结果 相当 。 而 温度 均 方 根 误差 在 近 地 层 达到 最 大 值 


1.88 KK， 其 中 在 800-1000 hPa 范围 小 于 ECMWF 再 分 
析 数 据 的 验证 结果 ， 而 在 其 余 压强 处 二 者 的 验证 精度 
相当 。 湿 度 平均 偏差 比 ECMWF 再 分 析 数 据 的 验证 结 


果 明 显 增 大 ， 


而 湿度 均 方 根 误 差 在 250 hPa 处 达到 最 


KIE 21.50%， 在 350-1000 hPa 范围 均 保持 在 16% 以 
内 ， 其 中 在 250~500 hPa 范围 内 的 湿度 均 方 根 误差 大 


于 ECMWF 再 分 析 数 据 的 验证 结果 。 


对 于 RAOB 


还 是 湿度 廓 线 ， 反 演 精度 比 其 他 两 个 数据 集 的 验证 结 
果 都 差 。 反 演 温度 廊 线 的 精度 在 接近 地 面 的 压强 范围 


数据 的 验证 结果 来 说 , 无 论 温度 廊 线 


内 最 差 , 为 3.43 K, 600-800 hPa 范围 内 的 精度 保持 在 


2~3 K 之 间 , 其 余 压 强 处 精度 保持 在 1~2 开 。 湿 度 廓 线 E v 5 
的 反 演 精度 在 850 hPa 附近 最 差 ， 为 25.48% 。 时 间 匹 100} f. 1 ee 
12 +> * 
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profilés 结果 来 说 ， 从 图 9%(b) 中 可 以 看 出 ， 反 演 后 的 湿度 平均 
4.3.4 = NCEP 预报 廊 线 的 精度 对 比 偏差 比 背 景 亡 线 的 平均 偏差 有 很 大 的 减 小 ， 最 大 处 可 


为 了 说 明 MWHTS 观测 资料 对 数值 天 气 预报 同化 WNA 15%， 在 湿度 探测 通道 起 主要 贡献 作用 的 
系统 的 重要 意义 ， 把 本 文 构 建 的 反 演 系统 的 反 演 结 250-750 hPa 的 压强 范围 ， 湿 度 反 演 结果 对 预报 廓 线 
与 NCEP 6 小 时 预报 廊 线 的 精度 进行 对 比 。 考 虑 到 的 精度 有 很 大 的 提高 ， 最 大 可 提高 18% 。 通 过 以 上 分 
NCEP GFS 数据 资料 的 存档 时 间 是 2015 年 1 月 , Ae 析 可 知 ， 本 文 针对 FY-3C/MWHTS 陆地 晴空 条 件 下 的 
利用 2015 年 1 月 到 3 月 的 晴空 NCEP 分 析 数 据 ， 取 观测 亮 温 所 建立 的 一 维 变 分 系统 ， 温 度 反 演 精度 对 预 
代 原 系统 中 的 ECMWF 再 分 析 数据 作为 先 验 信息 , 重 。” 报 廓 线 的 精度 改善 不 大 ， 而 湿度 反 演 结果 对 改善 预报 


新 生成 背景 协 方差 矩阵 B， 正 演 误差 协 方差 矩阵 E, 廊 线 精度 有 很 大 帮助 ， 从 而 说 明了 MWHTS 观测 资料 
以 及 逐 扫描 点 偏差 校正 方法 的 校正 系数 ， 生 成 方法 不 的 应 用 价值 以 及 对 数值 天 气 预 报 同 化 系统 有 重要 意 


变 。 而 背景 廊 线 仍然 选择 NCEP 6 “IN THRICE. DI 义 。 

NCEP 分 析 资 料 为 基准 验证 2015 年 4 月 到 5 月 的 ”4.3.5 与 AMSU 资料 的 反 演 结果 的 验证 分 析 
MWHTS 亮 温 反 演 结果 。MWHTS Aint NCEP 分 析 AMSU 的 通道 设置 的 功能 与 MWHTS 相似 , 同样 
数据 以 及 6 小 时 预报 数据 仍然 采用 0.5 小 时 的 时 间 误 有 对 大 气温 湿 廓 线 反 演 的 能 力 ， 其 观测 资料 已 业务 
差 ，0.05° 的 经 纬度 误差 进行 匹配 ， 按 照 4.1 节 的 反 演 化 应 用 。 选择 与 本 文 使 用 MWHTS 观测 资料 的 地 理 范 
WER, 共 得 到 102 组 反 演 结果 , 与 背景 廊 线 , BY NCEP 姑 和 时 间 范 围 都 相同 的 NOAA-18 Æ AMSU-A 和 
6 小 时 预报 廊 线 的 精度 对 比 ， 如 图 9 所 示 。 AMSU-B 陆地 观测 资料 按照 4.1 节 的 反 演 步 又 进行 反 


演 ， 通 过 反 演 结果 对 背景 廓 线 精度 的 改善 作用 ， 来 进 
一 步 评 估 MWHTS 观测 资料 。 其 中 ,为 了 同时 验证 背 
景 廊 线 对 有 反 演 精度 的 影响 ,选取 2014 Æ ECMWF ERA 
Interim 再 分 析 数 据 集 的 平均 廓 线 取 代 NCEP 6 小 时 预 
报 廓 线 作为 背景 廓 线 ， 针 对 MWHTS. AMSU-A fll B 
的 亮 温 数据 反 演 分 别 得 到 了 90、93 和 484 组 有 反 演 结果 ， 
利用 ECMWF 再 分 析 数 据 对 反 演 结果 进行 验证 ， 如 图 
10 所 示 。 
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图 10 与 AMSU 观测 资料 的 反 演 结果 对 比 


Fig. 10 Comparison of inversion results of MWHTS and AMSU 


与 图 7 相 比 , 从 图 10 中 可 以 看 出 , 当 使 用 气候 学 
数据 集 的 平均 廊 线 作为 背景 廊 线 时 ， 背 景 廓 线 的 均 方 
根 误 差 增 大 ， 导 致 了 温 湿 廓 线 的 反 演 精度 有 明显 的 下 
降 。 因 此 进一步 说 明了 在 本 文 构建 的 反 演 系统 中 ， 背 
景 廓 线 的 选择 对 反 演 精度 的 重要 影响 。 在 图 10 中 , 使 
] MWHTS 观测 资料 的 温度 廓 线 反 演 精 度 在 10~200 
hPa 范围 ， 对 背景 廓 线 的 精度 提高 幅度 明显 优 于 
AMSU-A 的 观测 资料 ， 而 在 其 余 压强 范围 ， 不 如 
AMSU-A 观测 资料 的 反 演 结果 对 背景 廓 线 的 改善 作 
。 对 于 湿度 反 演 的 均 方 根 误差 来 说 ，MWHTS 和 
AMSU-B 的 观测 资料 在 750~1000 hPa 和 450 hPa 附近 
范围 对 背景 廊 线 的 精度 改善 作用 相当 ， 而 在 其 余 压强 
范围 MWHTS 观测 资料 对 背景 廓 线 精度 的 提高 幅度 高 
UU QUE 

另外 , 分 别 选择 基于 AMSU-A 和 AMSU-B 陆地 
晴空 观测 资料 的 神经 网 络 反 演 和 物理 反 演 结果 与 本 文 
的 反 演 结 果 进 行 对 比 。 其 中 统计 反 演 方法 选择 神经 网 
络 反 演 (Karbou et al., 2005)， 物 理 反 演 方法 选择 MIRS 
反 演 系统 (Boukabara et al., 2011)， 这 两 个 方法 与 本 文 
的 最 大 区 别 在 于 都 对 地 表 进行 了 分 类 处 理 ， 其 中 这 两 
个 方法 所 使 用 的 AMSU-A 和 AMSU-B 的 观测 资料 分 
别 来 自 NOAA-15 和 NOAA-18。 反 演 结果 的 对 比 均 选 
择 与 RAOB 数据 的 验证 精度 , 温度 反 演 的 均 方 根 误差 
在 100 hPa、300 hPa、500 hPa、800 hPa 和 950 hPa 处 
进行 对 比 ， 湿 度 反 演 的 均 方 根 误差 在 300 hPa. 500 
hPa、800 hPa #11950 hPa 处 进行 对 比 ， 如 表 2 和 表 3 
所 示 。 
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表 2 温度 反 演 的 均 方 根 误差 对 比 


Table 2 Comparision of root mean square difference for the 


temperature 
压强 /hPa 100 300 500 800 950 
本 文 /K 178 1.47 1.96 2.59 3.09 


神经 网 络 K 1.82 241 408 410 3.68 
MIRS 系统 /区 240 180 180 280 3.10 


#3 湿度 反 演 的 均 方 根 误差 对 比 


Table 3 Comparision of root mean square difference for the 


humidity 
压强 /hPa 300 500 800 950 
本 文 /(%) 23.83 21.89 22.78 18.26 


神经 网 络 /(9%) 17.33 17.31 20.00 24.05 
MIRS 系统 /(%) 83.00 60.00 41.00 30.00 


从 表 2 中 可 以 看 出 ， 本 文 的 温度 反 演 精 度 在 高 空 
范围 明显 高 于 神经 网 络 方法 和 MIRS 反 演 系统 对 
AMSU 观测 亮 温 的 反 演 结果 ， 在 低空 范围 与 MIRS Js 
演 的 精度 相当 ， 而 高 于 神经 网 络 方法 的 反 演 精度 。 从 
表 3 中 可 以 看 出 ， 在 低空 范围 本 文 的 湿度 反 演 精度 好 
于 神经 网 络 方法 ， 而 在 高 空 范围 ， 不 如 神经 网 络 反 演 
方法 的 反 演 精度 。 然 而 MIRS 反 演 系统 的 湿度 反 演 精 
度 较 差 ， 但 是 该 反 演 系 统 在 云雨 天 气 条 件 下 对 温 湿 廓 
线 的 反 演 能 力 是 本 文 构建 的 反 演 系统 所 不 具备 的 。 通 
过 对 比 可 以 发 现 ， 虽 然 本 文 构建 的 反 演 系统 为 了 减少 
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计算 的 复杂 性 ， 在 计算 地 表 发 射 率 时 并 未 进行 详细 的 
分 类 处 理 , 但 是 针对 MWHTS 晴空 观测 资料 的 反 演 具 
有 较 高 精度 。 
5 结 论 


本 文 建立 了 基于 MWHTS 观测 亮 温 的 一 维 变 分 反 
演 系 统 对 晴空 条 件 下 的 陆地 上 空 大 气温 湿 廓 线 进行 反 
演 。 针 对 地 表 发 射 率 计算 复杂 性 高 且 精 度 较 差 的 问题 ， 
本 文 从 先 验 信息 的 角度 出 发 ， 在 反 演 过 程 中 首次 提出 
了 背景 协 方差 矩阵 的 联合 矩阵 和 单独 矩阵 对 温 湿 廓 线 
进行 组 合 反 演 的 方法 ， 有 效 的 减 小 了 反 演 过 程 的 复杂 
性 ， 且 反 演 性 能 优 于 根据 地 表 类 型 进行 分 类 计算 的 反 
演 方法 。 对 反 演 结果 利用 多 个 数据 源 进 行 反 演 精度 的 
验证 ， 证 明了 结果 的 可 靠 性 。 
利用 本 文 所 构建 的 反 演 系统 分 别 对 MWHTS 和 
AMSU 的 亮 温 数据 进行 反 演 ， 对 于 温度 反 演 来 说 ， 
MWHTS 的 反 演 精度 在 10~200 hPa 的 高 空 范 围 优 于 
AMSU 的 反 演 结果 ， 而 对 于 湿度 反 演 来 说 , MWHT 
在 除了 750~1000 hPa 和 450 hPa 附近 范围 与 AMSU 有 具 
有 相当 的 反 演 精度 外 , 其 他 压强 范围 均 优 于 AMSU 的 
反 演 性 能 。 反 演 结果 的 对 比 表 明 MWHTS 在 高 空 范围 
有 
资 


Nn 


有 有 较 强 的 温度 探测 能 力 ， 而 在 整个 压强 范围 具 
较 强 的 水 汽 探测 能 力 。 从 而 也 证 明了 MWHTS 观测 
料 具 有 较 高 的 质量 。 同 时 ，MWHTS 观测 资料 的 反 演 
结果 对 NCEP 6 小 时 预报 廓 线 精度 的 改善 作用 表明 了 
MWHTS 观测 资料 对 数值 天 气 预报 同化 系统 具有 重要 
的 意义 。 

然而 ， 在 与 各 个 数据 源 进 行 精度 验证 时 发 现 ， 背 
景 廓 线 的 误差 与 反 演 误差 的 分 布 一 致 性 较 高 ， 是 影响 
反 演 精度 的 重要 误差 来 源 ， 为 了 进一步 提高 系统 的 反 
演 精 度 ， 寻 找 高 精度 的 背景 廊 线 数据 源 是 下 一 步 的 工 
作 方 向 ， 同 时 模拟 亮 温 与 观测 亮 温 之 间 的 偏差 校正 效 
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果 对 反 演 精度 具有 重要 的 影响 ， 因 此 偏差 校正 方法 的 
优化 也 是 下 一 步 的 工作 重点 。 
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Obiective: The land surface emissivity has great effect on the atmospheric temperature and humidity sounding from 


space. Since the complication of land surface, the calculated emissivity is generally with low accuracy. In order to 


improve retrieval accuracy and reduce the computation complexity in retrieving atmospheric temperature and 


humidity profiles over land in clear sky, an one-dimensional variational retrieval system was built using measurements 
of microwave humidity and temperature sounder onboard the Fengyun-3C satellite (FY-3C/MWHTS). 

Method: By analyzing the priori information affecting the accuracy of inversion, A hybrid retrieval approach based on 
united matrix and individual matrix of background covariance is proposed. The method established a better 
correlation relationship between temperature and humidity profiles, reduce the error propagation of the retrieval 
temperature and humidity, and avoid the complicated land emissivity calculations according to different surface types. 
In addition, based on the correlation between the observed values from FY-3C/MWHTS and those simulated by the 
forward radiative transfer model, a statistical regression method is adapt in pixel by pixel correction procedure to 
correct the bias between the observed values and those simulated values. 

Result: This retrieval system retrieves temperature and humidity profiles over a part of China’ land in clear sky and 
validates the retrieval results with respect to the European Centre for Medium Range Weather Rorecasts(ECMWF) 
reanalysis data, National Centers for Environmental Prediction(NCEP) analysis data and  Radiosonde 
Observation(RAOB) data, respectively. With respect to the ECMWF reanalyse data, the maximum root mean square 
errors of the resulted temperature and relative humidity are 2.59K and 11.87% respectively. With respect to the NCEP 
analyse data, the maximum root mean square errors of the resulted temperature and relative humidity are 1.88K and 
21.5096, respectively. With respect to the RAOB data, the maximum root mean square errors of the resulted 
temperature and relative humidity are 3.43K and 25.48%, respectively. Comparison of the retrieval results with those 
of measurements of AMSU using the physical retrieval method and statistical retrieval method show that the MWHTS 
has higher accuracy. Comparison of the retrieval results with NCEP 6 hour forecast profiles show that the retrieval 
humidity profiles can improve the accuracy of forecast profiles, the retrieval temperature profiles can improve the 
accuracy of forecast profiles in the upper atmosphere. 

Conclusion: The proposed hybird approach using united matrix and individual matrix of background covariance can 
get satisfactory retrieval accuracy, although the land emissivity is calculated without classing the surface types in the 
retrieval system built in this paper. The retrieval system whose retrieval results are evaluated by root mean square 
errors with respect to three data sources has high accuracy and high reliability. Comparison of the retrieval results with 
those of measurements of AMSU suggests that MWHTS has a greater ability to probe the temperature in the upper 
atmosphere and humidity in the entire atmosphere. Comparison of the retrieval results with NCEP 6 hour forecast 
profiles suggests that MWHTS can get high-quality data for sounding atmosphere and is great importance to the 
numerical weather prediction radiance assimilation. 

Keywords: microwave humidity and temperature sounder(MWHTS), one-dimensional variational retrieval, 
temperature and humidity profiles, pixel by pixel correction, forecast profiles 
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